Limnology and Oceanography: Methods

The morphology of cyanobacterial colonies plays a key role in harmful cyanobacterial blooms, with implications for their vertical migration, resistance against grazing, and light availability. In this study, we introduce the use of Optical Coherence Tomography (OCT) to investigate the three-dimensional morphology of cyanobacterial colonies. The technique enables non-invasive 3D imaging of colonies up to several millimeters in size, providing access to detailed mesoscale morphological features. Gas vesicles inside cells were shown to strongly improve image quality. We describe the sample preparation and image acquisition protocol, as well as an image processing pipeline that extracts mesoscale morphological features and provides a volumetric visualization of colonies. The method was tested for representative colonies of different cyanobacterial species while a dataset of volumetric images and measured mesoscale features was acquired for natural colonies of Microcystis. We demonstrate the utility of 3D imaging by quantifying the effects of irregular colony morphologies on their flotation velocity and the light availability within colonies. We anticipate OCT to become a key imaging technique to monitor populations of cyanobacterial colonies and investigate colony formation, with potential extensions to other colonial and aggregated organisms in freshwater and marine environments.

Passive sampling is the deployment of a collection material in the environment to continuously capture environmental DNA (eDNA) over time, offering the potential to integrate biodiversity signals while reducing the need for repeated active water collection. However, the mechanisms governing eDNA capture and retention on passive samplers remain poorly understood, limiting the interpretation of passive eDNA signals and their broader application. Here, we investigated the mechanistic performance of glass fibre passive samplers using controlled mesocosm experiments with three invasive freshwater bivalves: zebra mussels (Dreissena polymorpha), quagga mussels (Dreissena bugensis), and Asian clams (Corbicula fluminea). Specifically, we quantified eDNA accumulation dynamics, evaluated the contribution of different eDNA states, tested the persistence of captured eDNA, and compared passive sampler signals with conventional active sampling. Passive samplers rapidly accumulated target eDNA within hours of deployment, after which concentrations either plateaued or continued to increase depending on species. Sequential transfer of passive samplers between mesocosms containing different species showed that previously captured eDNA declined while new target eDNA accumulated to concentrations comparable to freshly deployed samplers, demonstrating continual turnover rather than permanent retention. Dissolved eDNA showed little evidence of accumulation beyond the concentration retained in the pore water within the membrane, suggesting that it is unlikely to be the dominant contributor to long-term passive sampler signals. Instead, the observed variability among replicate samplers, together with the physical properties of glass fibre membranes, suggests that membrane-bound and particulate eDNA are the primary contributors to passive eDNA capture. Collectively, these findings support a model in which glass fibre passive sampler signals reflect a dynamic equilibrium between ongoing eDNA capture and concurrent loss processes rather than cumulative accumulation over time. This mechanistic framework provides a foundation for interpreting passive eDNA data and informs the future development of passive sampling materials, deployment strategies, and biodiversity monitoring applications.

The picocyanobacterium Prochlorococcus is a fundamental contributor to ocean primary productivity. While its free-living population has been extensively studied, primarily using flow cytometric analyses, the size and distribution of its particle-associated population is not well understood. Using filter fractionated samples from cruises in the Pacific Ocean, Atlantic Ocean, and Mediterranean Sea, we generated metagenomic data using internal standards, yielding absolute genome equivalent counts of Prochlorococcus cells in different size fractions. We used these data to model a relationship between relative and absolute genome equivalent counts, yielding a correction factor that we validated using published datasets. We then applied the correction factor to size-fractionated global metagenomic data from the TARA Oceans Project, which has widespread Prochlorococcus cells in size fractions >1.6 m throughout the transects, to calculate the fraction of the total Prochlorococcus population in large size fractions. The ''particle-associated'' population fraction increased with net primary productivity. Dissolved inorganic carbon was also directly correlated with increased particle association, which, combined with other evidence, could indicate an association with upwelling. We also examined the relationship between particle-associated population and carbon export at 150 m by incorporating published estimates of carbon flux based on TARA optical scattering data. This study highlights the potential importance of particle-associated Prochlorococcus to carbon flux in marine ecosystems and offers a way to convert relative to absolute genome equivalents of microorganisms in archival metagenomic datasets.

Diatoms are major contributors to marine primary production, yet current approaches for monitoring their abundance and function rely on coarse satellite chlorophyll estimates or sparse cell count and carbon fixation measurements. Molecular markers are a promising approach for high-resolution measurement of both abundance and metabolic activity through analysis of environmental DNA (eDNA) and RNA (eRNA). We present an isothermal quantitative recombinase polymerase amplification (qRPA) assay targeting rbcL gene copies and transcripts of marine diatoms, operating at low temperature and producing results in less than 15 min. We demonstrate specificity and calibration across diverse diatom taxa, then apply the assay to eDNA and eRNA samples from the Mare Chiara Long-Term Ecological Research site in the Bay of Naples, Italy, alongside microscopy, chlorophyll, physicochemical, and carbon-fixation data. Diatom rbcL DNA tracked abundance across five orders of magnitude despite seasonal shifts in community composition. Combining molecular and optical data revealed increased cellular rbcL copies and chlorophyll in low-light winter populations, suggesting enhanced photosynthetic capacity despite lower abundance. Furthermore, rbcL RNA reflected total carbon fixation rates and identified populations with differing carbon fixation activity. These results support rapid, RPA-based rbcL quantification as a robust approach for biomolecular ocean observing.

Ranges of small benthic fauna are notoriously difficult to assess. In some of these cases, modern eDNA methods can shed light on species occurrence. Here we conduct an exploratory study on the fish eDNA recoverable from the gastrovascular cavities of the easy-to-sample pore water siphoning benthic invertebrate, Cassiopea, across six sites within the Florida Keys. Twenty-seven fish 12S identities were recovered from water samples, two from sediment samples, and seventeen from Cassiopea gut swabs. In total, thirty-two different species were identified from nineteen families, including one shark species (Ginglymostoma cirratum), and five species of cryptobenthic reef fishes (f: Gobiidae, Labrisomidae). Additionally, five species were identified from medusae samples that were not recovered in water or sediment samples. The species identities recovered may provide insight into the fish in direct proximity to Cassiopea assemblages, as well as indicate that Cassiopea may accrue disproportionate eDNA from cryptobenthic reef fish compared to surrounding environmental samples. The unorthodox sampling technique of using eDNA recovered from jellyfish stomachs yields another avenue for epibenthic community data acquisition.

Coral reef monitoring needs scalable, non-invasive tools to complement resource-intensive traditional survey methods. Passive Acoustic Monitoring (PAM) offers a promising supplement, but its effectiveness is limited by the difficulty of attributing recorded sounds to species outside of previously well-characterized taxa. Using Omnidirectional Underwater Passive Acoustic Cameras (UPAC-360), we identified sounds from 31 reef fish species across 14 families on the Kona coast of Hawaii Island, including 13 not previously documented as soniferous. By releasing video and audio specimens, we have created the largest open-access collection of in-situ reef fish sounds to date for the Pacific. A subset of acoustically distinctive taxa--such as Hawaiian Dascyllus (Dascyllus albisella), Lei Triggerfish (Sufflamen bursa), soldierfishes (Myripristis spp.), wrasses, and herbivorous grazers--were identifiable in PAM recordings through manual acoustic and spectrogram review. Through identifying particular sounds linked to species with different ecological roles, these sounds have the potential to serve as indicators of reef function to increase the information and value coming from PAM surveys of Hawaiian and Pacific coral reefs.

O_LIEnvironmental DNA (eDNA) methods are developing rapidly for ecological surveys, and passive eDNA sampling has emerged as a promising approach for integrating DNA signals over deployment time. However, how deployment duration affects the amount of detectable DNA retained by a sampler remains poorly understood. C_LIO_LIHere, an analytical model was developed to examine how DNA input, degradation, finite substrate capacity and residual retention of degraded DNA shape passive eDNA accumulation. The model distinguishes detectable adsorbed DNA from degraded, non-detectable DNA that may remain on the substrate and continue to occupy capacity. The residual-retention parameter,{theta} , represents the fraction of degraded DNA that remains capacity-occupying, with{theta} = 0 corresponding to complete replacement and{theta} = 1 to complete non-replacement. C_LIO_LIThe model predicts three key behaviours. First, when degraded DNA does not occupy substrate capacity ({theta} = 0), detectable eDNA accumulates monotonically towards equilibrium, but equilibrium recovery increases less than proportionally with DNA input. Thus, passive-sampler measurements can compress quantitative differences in environmental DNA supply. Second, when degraded DNA remains capacity-occupying ({theta} > 0), detectable eDNA can reach a finite peak and subsequently decline. Higher DNA input increases peak yield but shifts the peak earlier, whereas greater substrate capacity increases peak yield and delays the peak. Third, under prolonged deployment with{theta} > 0, a higher-input condition can yield less detectable eDNA than a lower-input condition, reversing the expected input-rate ranking. C_LIO_LIThese results show that passive eDNA recovery can follow saturating, unimodal or intermediate dynamics depending on substrate capacity and post-adsorption DNA fate. Thus, retrieval time cannot be optimised by adjusting deployment duration alone. Although investigators can choose deployment duration and sampler design, including substrate capacity, optimisation also requires calibration or explicit assumptions about ambient DNA supply, DNA degradation rate and residual retention of degraded DNA. C_LI

Coral reef fish assemblages are widely used as indicators of ecosystem condition, yet manual annotation of underwater video remains a major bottleneck for scalable biodiversity monitoring. Despite rapid progress in automated detection, ecologically realistic and publicly available datasets remain scarce, particularly for the Western Indian Ocean. Here, we present WIO-ReefFish, a reef fish detection dataset derived from diver-operated line-intercept transects and designed for ecological monitoring under natural survey conditions. WIO-ReefFish comprises 1,000 ultra-high-definition images (3840 $\times$ 2160 pixels) and 6,768 exhaustive bounding-box annotations spanning 24 taxonomic categories, thereby preserving full-frame assemblage structure in complex reef scenes. We also establish a standardized benchmark across nine object detection models under two complementary protocols: class-aware detection and class-agnostic fish localization. Detection performance was consistently higher under the class-agnostic protocol. The best-performing model (RT-DETR) improved from 0.48 mAP50 in the class-aware setting to 0.70 mAP50 when taxonomic constraints were removed, indicating that taxonomic discrimination remains substantially more challenging than fish localisation in reef imagery. Spatially independent evaluation revealed a pronounced generalisation gap, particularly for taxonomic detection, whereas class-agnostic fish localisation remained substantially more robust across transects and countries. Together, these results establish WIO-ReefFish as a realistic benchmark for automated reef fish detection and provide a foundation for more robust computer-vision tools in coral reef biodiversity monitoring. The WIO-ReefFish dataset and associated benchmarking resources are publicly available.

This brief communication presents two solutions for calculating absolute measures of error from time-difference-of-arrival (TDOA) positioning in underwater acoustic telemetry arrays. First, a Monte Carlo estimation of TDOA positioning error is derived. Next, a computationally inexpensive, approximate solution to the Monte Carlo method is presented. This approximate solution is achieved by solving the Jacobian of a closed-form TDOA positioning model. The positioning error covariance matrix returned from either method can then be used to report the accuracy of TDOA positions or utilized in state-space positioning models. Finally, calculations of the expected radial error are shown which serves as a simple summary statistic for reporting positioning error in real units.

Offshore wind farms are rapidly expanding to meet growing demands for renewable energy, with development expected to extend further offshore into deeper waters. This expansion requires a robust understanding of the long-term ecological consequences of offshore wind farms (OWFs) and how these may interact with ongoing climate change. We used the coupled hydrodynamic-ecosystem-biogeochemical water-column model (GOTM-ERSEM-BFM) to investigate ecosystem-wide responses to environmental changes associated with OWFs and climate warming. Specifically, we examined OWF-related scenarios of reduced benthic suspension-feeding activity, representing potential effects of contaminant emissions from OWFs, and reduced wind forcing, together with increased sea surface temperature. The scenarios were simulated individually and in combination to explore potential interactive effects. These scenarios were simulated at two contrasting locations in the North Sea, representing a well-mixed coastal site and a seasonally stratified offshore site. The coastal site exhibited comparatively modest ecosystem responses across the scenarios, whereas responses were generally stronger at the deeper offshore site. At the offshore site, changes in stratification altered vertical nutrient dynamics and contributed to pronounced differences in ecosystem responses between the surface and bottom layers. Our results demonstrate that ecosystem responses to OWF-related and climate-driven environmental changes are strongly dependent on local environmental conditions, suggesting that ecological consequences may differ substantially as wind farm development expands into deeper offshore environments.

Marine diatoms contribute to about 20% of global primary production. We present the first basin-scale, multiyear assessment of diatom communities in the North Pacific, combining taxonomically high-resolution RuBisCO large subunit gene (rbcL) metabarcoding with concurrent environmental measurements. Using a nine-year time series of daily samples resolved at the species level via ~500 bp rbcL fragments, we performed multivariate analyses across biogeographic provinces, identifying significant correlations between community structure and environmental drivers such as temperature and macronutrient availability. We report the prevalence of a previously overlooked centric diatom species in the North Pacific, Eunotogramma lunatum, which appears to be near-dominant even in subarctic high-nitrate, low-chlorophyll waters where pennate diatoms are typically favored. These results demonstrate the power of rbcL for large-scale ocean monitoring and provide a critical baseline for future studies of diatom population dynamics, climate change impacts, and ecosystem resilience in a key marine region.

Host genetic markers often fail to resolve regional origins in highly connected or panmictic marine species. The Japanese jack mackerel, Trachurus japonicus, is a commercially important fishery species around Japan that shows little or no detectable population structure. Here, we used nationwide multi-omics profiling to compare host genomic variation and gut microbiome composition in wild T. japonicus collected from coastal regions across Japan. We generated MIG-seq data for 43 individuals and 16S rRNA gene profiles for 24 individuals; after quality filtering, 19 individuals remained for matched host-microbiome comparison. Genome-wide host SNP analyses showed weak or absent geographic population structure, consistent with previous evidence of panmixia in Japanese waters. In contrast, gut microbiome composition showed geographic structuring based on Bray-Curtis dissimilarity and PERMANOVA, and this pattern was not explained by proximity to river mouths or host-related variables. Locality- or individual-associated bacterial lineages contributed to the observed differences in the microbiome, while chloroplast-associated and Cyanobacteria-assigned ASVs suggested recent dietary or environmental input. These results indicate that gut microbiome can show regional biological variation not apparent from host genetic markers alone. Our study provides a proof-of-concept example of integrating host genomics and gut microbiome profiling to evaluate regional characteristics and origins in highly connected marine animals.

This research addresses a central question in Palaeolithic research: how hunter-gatherer mobility was structured across space and time in the Cantabrian Region (northern Iberia), which has human occupation evidence spanning from the Middle Pleistocene through the Holocene. A multidisciplinary framework integrating primarily {delta}3S isotope values, combined with {delta}{superscript 1}3C and {delta}{superscript 1}N, palaeoproteomics, Bayesian age modelling, palaeoclimatic reconstruction, isoscape mapping, and ecological diversity was developed. A total of 905 animal bone collagen samples, with evidence of anthropogenic modifications, from 16 key archaeological sites from the Mousterian to Mesolithic (Marine Isotopic Stage 5 to 1, between 100-7 ka BP) were analysed, permitting the reconstruction of spatial patterns of resource exploitation and human mobility. The {delta}3S isotope values show weak, inconsistent relationships with climatic proxies, suggesting that sulfur signatures are primarily driven by geographic and ecological factors rather than climate. Strong spatial trends are observed, with higher {delta}3S values in coastal zones and lower values inland. Diachronic trends reveal marked shifts in human mobility: smaller ranges during the Mousterian, increasing mobility through the Chatelperronian and especially the Aurignacian, followed by reduced mobility in the Gravettian and Solutrean, and renewed territorial expansion during the Magdalenian and, likely, the Azilian. In contrast, the Mesolithic is characterised by decreased mobility and thus increased territoriality in both coastal and inland contexts. Faunal isotope values and isoscape predictions reveal that some animals were acquired beyond local foraging ranges during the Palaeolithic, particularly in inland regions with lower {delta}3S values. Isotopic niche analyses indicate partial interspecific overlap consistent with ecological flexibility. Macromammal and micromammal diversity exhibit contrasting patterns, with a significant negative correlation in Simpson and Shannon indices. Macromammal diversity correlates negatively with {delta}3S values, linking increased hunting diversity to expanded catchment areas and longer-distance foraging, whereas micromammal diversity shows positive correlations with {delta}3S, {delta}{superscript 1}3C and {delta}{superscript 1}N reflecting stronger climatic influence. Overall, these results demonstrate that hunter-gatherer behaviour in northern Iberia during the Middle and Late Palaeolithic was highly dynamic, combining logistical and residential strategies that shifted in response to changing environmental conditions, resource distributions and cultural adaptations.

The marine biodiversity inhabiting rocky shores in the Colombian Pacific remains largely undocumented, primarily due to geographic isolation, logistical challenges, and socio-political constraints. To address the existing knowledge gap, we conducted an expedition to enhance baseline biodiversity knowledge in rocky shores by integrating multiple complementary approaches, including visual censuses, specimen collection with morphological identification, environmental DNA (eDNA) metabarcoding and DNA barcodes. eDNA samples were collected at four coastal sites adjacent to rocky substrates, along with biological specimens obtained from fourteen locations through SCUBA diving at depths ranging from 1 to 25 meters. Tissue samples were subjected to genomic DNA isolation, followed by the generation and validation of cytochrome c oxidase subunit I (COI) barcode sequences, which were subsequently corroborated through taxonomic assessment to ensure accurate species identification. eDNA metabarcoding analyses yielded over 7 million high-quality sequence reads. Although taxonomic resolution at the species level was constrained by the limited completeness of reference sequence databases, a total of 106 species and 83 families were successfully identified, predominantly within the classes Actinopteri, Chondrichthyes, and marine mammals. From the 769 specimens obtained we generated 871 sequences, including 414 validated COI barcodes representing 76 species across 64 families. The integration of DNA barcoding and eDNA approaches resulted in over 1,400 taxonomic detections spanning five phyla, with only six species shared between methodologies. Richness and diversity varied among sites, and revealed significant differences along the coastline between Jurado and Cupica Gulf. All sequences were deposited in BOLDsystems database under the CCBIO project and were visualized through OBIS and GBIF databases. These findings provide the first molecular-based baseline for rocky shore biodiversity in the Colombian Pacific, highlighting the value of integrative approaches for monitoring and conservation.

Deep learning tools are increasingly used today, particularly in medical segmentation. A gap nonetheless remains in automating segmentation for insects. This work addresses the following question: can a generalist segmentation model, trained on several phylogenetically related orthopteran species, reliably automate head tissue segmentation from micro-CT images? To answer this, we used nnU-Net, a self-configuring 3D deep learning segmentation framework originally developed for medical imaging, whose core function, learning to recognize tissues of interest, applies directly to this context. Six anatomical classes were automated, comparing two training strategies: sequential fine-tuning, which adds species one at a time under the assumption that progressive learning would strengthen predictive power, and from-scratch training, in which the model learns the entire dataset simultaneously. The fine-tuning model (ModelB) reached a Dice coefficient (a measure of overlap between automated segmentation and manual ground truth, ranging from 0 to 1) of 0.7715, compared to 0.7664 for the from-scratch model (ModelC). Although both models produced accurate automated segmentations, no significant difference was found between the two training strategies (paired Wilcoxon test, n = 24, p = 0.243). Despite a dataset limited to 20 individuals and the absence of one method clearly outperforming the other, the models remain usable across the three species studied (Gryllus bimaculatus, Loxoblemmus equestris, L. doenitzi), including in the presence of pronounced sexual dimorphism. It reduces a 20 hour segmentation task to under a minute.

Dense semilandmarks describe 3D surfaces with hundreds to thousands of points, and sliding them by bending energy or Procrustes distance is a near-universal default. Three questions remain open: does dense sampling add shape beyond fixed landmarks, how many points are needed, and does sliding help or harm? Real specimens cannot answer them: the true correspondence is unknown. We tested two workflows, ALPACA (single-template registration) and DeCAL (landmark-anchored correspondence), on 496 mouse skulls at 250-1,000 points, with and without sliding, scored by surface reconstruction. We repeated it on 500 synthetic skulls with exact correspondence, measuring each point's distance to its true homologue. Dense semilandmarks lowered error for almost every specimen; the fixed landmarks added little but supplied anchoring the semilandmarks could not, and the anchored method was more accurate. The benefit saturated near 250 points for ALPACA but kept improving to 1,000 for DeCAL. Procrustes-distance sliding harmed every configuration; bending-energy sliding helped only a poor, landmark-free correspondence, vanishing once anatomical anchors spanned the form. Match the sliding decision to the correspondence in hand: relax a poor one, leave a good one alone, never slide toward the mean. Known-correspondence specimens offer a general test of landmarking and sliding against ground truth.

Full-length 16S rRNA gene sequencing using Oxford Nanopore Technologies has emerged as a promising approach to improve species-level resolution in microbiota studies. However, the accuracy of taxonomic assignment remains highly dependent on the bioinformatics s used to process Nanopore long-read data. Therefore, the only way to ensure a good level of certainty in obtained results is to use positive controls in the form of mock communities in the experimental designs. In this study, we compared the performance of Epi2Me 16S (using Minimap2 or Kraken2) workflows provided by Oxford Nanopore Technologies and an EMU workflow for full-length 16S rRNA gene analysis. Using a commercial mock community sequenced across multiple Nanopore runs, taxonomic assignment accuracy and reproducibility was evaluated. Epi2Me-Kraken2 exhibited 18 % of incorrect genus-level assignments and failed to identify 3 species present in the mock community. While Epi2Me-Minimap2 achieved an excellent genus-level classification, reporting 9 % of sequences assigned to a genus not in the mock community, species-level assignments were inconsistent for several community members such as Listeria. In contrast, EMU provided accurate and consistent species-level taxonomic profiles, with all species correctly identified while keeping the number of genus absent from the mock community at 1.2%. ImportanceThese results highlight that Epi2Me integrated workflows are not the best option for specie-level taxonomic assignation. More importantly, this paper underscores the importance of routine inclusion of positive controls for microbiota studies, in the form of mock communities, as a critical safeguard for accurate data interpretation. Without the use of a mock community, a paper published would be at risk of reporting wrong observations and inaccurate conclusions.

Sexual coral propagation is an emerging technique capable of producing large numbers of corals for coral transplantation and reef rehabilitation. In contrast to asexual coral propagation, sexual propagation increases genotypic diversity and can be used for selective breeding to enhance coral heat tolerance or other desirable traits. However, implementation at meaningful ecological scales is hindered by high mortality during early life stages, high costs associated with nursery rearing facilities, and labour-intensive outplanting methods. To overcome these issues, we developed the CoralAssist Plug (CAP), a ceramic device designed for the rapid and cost-effective outplanting of sexually propagated corals in large numbers that maximises post-outplant survivorship. CAPs combine three important functional features: 1) built-in microrefugia to protect juvenile corals from grazing, 2) a relatively small size, 3 by 1 cm, that is easy to handle and stack efficiently without compromising the survivorship of corals, and 3) a hole in the middle that facilitates handling and attachment. CAPs were settled with Acropora aff. digitifera and outplanted to a reef crest after 1 to 6 months of ex situ nursery rearing. A 3-person dive team was able to outplant ~120 CAPs in one 90-minute shallow dive (just over 2 minutes per CAP per person). With longer nursery durations of 6 months, it was possible to achieve 36 % yield (i.e., the proportion of devices with a surviving coral) 4-years post-outplant. With nursery durations shortened to 1 month, we were able to attain 24 % yield 3-years post-outplant. Microrefugia significantly enhanced post-outplant survivorship leading to an 11 % increase in yield 4 years post outplant compared to devices without microrefugia. Outplanted corals that had reached adult size, were self-attached and were reproductively mature after 4 years. Our results suggest that CAPs can play a meaningful role in reef rehabilitation by efficiently introducing sexually propagated corals into natural populations with clear applications to assisted evolution techniques, such as selective breeding.

Antarctica represents a significant, unresolved, and unstable source of methane to the atmosphere. To advance our understanding of the biological filter of methane in Antarctica, here we identify the taxa and functional genes present during methane oxidation in an Antarctic Methane Seep. Methane oxidation was present in all sediments, including in a non seep control site. Using 16S rRNA analysis alongside metagenomics, we found that ANaerobic MEthane oxidizing (ANME) archaea coupled to Sulfate-Reducing Bacteria (SRB), documented as the most important marine methane sink in other locations, were not present. Instead, we observed the presence of denitrification-dependent methane oxidizers, including the anaerobic genus Candidatus Methylomirabilis, alongside the nitrate reducing archaea Candidatus Methanoperedens through short-read metagenomic classification. In addition, we note the presence of multiple aerobic methanotrophs, with a particularly high abundance of the Methylobacter, Methylomonas, and Methyloprofundus genera. Our results support denitrification-mediated methane oxidation and aerobic methanotrophy as the primary potential methane sinks in the Ross Sea. The widespread methane oxidation, including in control sediment, combined with the possibility of anaerobic methane oxidation linked to denitrification rather than sulfate reduction highlights the ubiquity and uniqueness of the Antarctic methane cycle.

Freshwater bioassessment relies on assessing aquatic assemblages to infer ecological conditions, yet conventional surveys require extensive field sampling, specimen processing, and specialized taxonomic expertise. Existing environmental DNA (eDNA) methods have not yet provided a practical alternative to conventional macroinvertebrate assays in part because current approaches cannot feasibly recover broad taxonomic diversity at sufficient taxonomic resolution. Here, we evaluated targeted hybridization capture of mitochondrial cytochrome oxidase I (COI) target sequences as a unified molecular approach for cross-phylum freshwater bioassessment. Environmental DNA was collected at 18 sites along 63 km of Boulder Creek spanning nearly 1,500 m of elevation from forested headwaters to agricultural plains. COI targets were enriched using custom RNA bait panels designed to target regional freshwater arthropods, annelids, and molluscs. Hybridization capture increased recovery of COI sequences [~]1,760-fold relative to unenriched shotgun libraries, generating Folmer-region COI contigs that averaged [~]400 bp. Across the watershed, we recovered sequences for approximately 450 macroinvertebrate genera across 8 phyla. Detected macroinvertebrate richness averaged 56 genera per site and increased down Boulder Canyon before declining downstream of the city. Macroinvertebrate assemblage composition from hybridization capture paralleled patterns observed with past conventional bioassessment. These results demonstrate that targeted hybridization capture enables robust, cross-phylum detection of species used for freshwater bioassessment from environmental DNA.